Ammonia volatization from conventional and stabilized fertilizers, agronomic aspects and microbiological attributes in a Brazilian coffee crop system

We aimed to quantify the N losses through volatilization of the main conventional and stabilized N fertilizers applied in coffee plantations. Additionally, we also assessed microbiological attributes of the soil (microbial biomass carbon (MBC); microbial biomass nitrogen (MBN); microbial basal respiration (MBR); metabolic quotient (qCO2); urease, β-glucosidase, acid phosphatase, and arylsulfatase activities) and agronomic aspects of the crop (N content in the leaves and beans, yield, and N exportation by the beans). Treatments consisted of the combination of three fertilizers (ammonium nitrate - AN, conventional urea - U, and urea with N- (n-butyl) thiophosphoric triamide (NBPT) - UNBPT, and five doses of N (0, 150, 275, 400, and 525 kg ha-1 year-1 of N), with four replicates, totalling 60 experimental plots. In the two crop seasons evaluated, daily and cumulative losses of N-NH3 from the split fertilizer applications were influenced by the N fertilizer technologies. The application of U resulted in losses of 22.0% and 22.8% for the doses of 150 and 400 kg ha-1 year-1 of N. This means that 66 and 182 kg ha-1 of N-NH3 were lost, respectively, at the end of six fertilizations with U. UNBPT reduced urease activity and N-NH3 losses compared to conventional urea, avoiding the volatilization of 15.9 and 24.3 kg ha-1 of N. As for AN, N-NH3 losses did not exceed 1% of the applied dose, regardless of the weather conditions during the fertilization. Urease activity was higher on days of maximum NH3 volatilization. There was an effect of the N sources (NS), soil sampling time (ST), and their interaction (NS × ST) on the MBN and arylsulfatase activity. The N sources also influenced the MBC and the qCO2. A substantial amount of N was removed from the system by the beans and husks of the harvested fruits. Our study showed that N fertilizer technologies are interesting options to reduce N-NH3 losses by volatilization, increase N retention in the soil, and improve microbiological attributes and the sustainability of coffee production systems.


Introduction
Brazil is the world's largest producer and exporter of coffee (Coffea arabica L.).Increasing fertilizer use efficiency through innovative technologies reduces costs of production, improves the quality of the beverage, and makes farming activity more sustainable and competitive, especially in times of high input prices.The application of nitrogen (N) fertilizers is essential for the productivity of this high-value crop.Brazilian agriculture consumed 9.2 million tonnes of N fertilizers in 2020 (FAOSTAT, 2021), in which coffee production alone was estimated do consume 10-15% of the country's N (1.4 million tonnes per year).
The intense dynamics of the N transformations in the soil and the diverse pathways it can be lost in coffee production systems result in low N fertilizer use efficiency (NUE) (Salamanca-Jimenez et al., 2017).Studies using 15 N demonstrated that coffee plants recover less than 25% of the N fertilizer when conventional urea is used (Cannavo et al., 2013;Salamanca-Jimenez et al., 2017), a complex scenario considering that approximately 50% of the global N fertilizer production is represented by urea (Heffer and Prud'homme, 2016;FAOSTAT, 2021).
The enhanced-efficiency fertilizers (EEFs) encompass many technologies for fertilizers and its one of the most studied strategies to increase NUE (Chien et al., 2009;Azeem et al., 2014;Timilsena et al., 2015) and mitigate greenhouse gas emissions in agriculture.For N fertilizers, these technologies are divided in the following main categories: conventional fertilizers, stabilized fertilizers, slow-release fertilizers, controlled-release fertilizers, and their blends (Guelfi, 2017).The proper use of EEFs is an interesting option to reduce N losses, carbon footprint, and climate change, in addition to increasing N retention in the soil (Snyder, 2017;Zhang et al., 2021;Freitas et al., 2022).
Stabilized fertilizers have additives, such as urease inhibitors (NBPT, NPPT, 2-NPT, and Duromide) and nitrification inhibitors (DCD, DMPP, and DMPSA), that inhibit or delay a determined stage of the N transformation process in the soil (Guelfi, 2017;Byrne et al., 2020).In Brazilian coffee cultivation, urease inhibitors are already widely used, while nitrification inhibitors or the combination of these two subgroups are still not significant in the national market.Currently, some authors have highlighted the need to evaluate the possible negative impacts of some additives on human health in future studies (Cai and Akiyama, 2017), which may be related to the persistence of some molecules in the environment.
The N-(n-butyl)thiophosphoric triamide (NBPT) is currently the most widely used molecule as a urease inhibitor (Sanz-Cobena et al., 2011).NBPT is reported to effectively reduce N-NH 3 losses through volatilization when added to urea (Adotey et al., 2017;Santos et al., 2020).The reduction of NH 3 losses is related to the maintenance of the N in the form of amide through a temporary inhibition of the urease activity, the enzyme which catalyzes the hydrolysis of urea (Santos et al., 2021;Freitas et al., 2022).This is due to the ability of the NBPT to oxidize into its analog, NBPTO, which forms stable complexes with the enzyme, causing its inactivation (Sanz-Cobena et al., 2011).The potential inhibition of the additive depends on the concentration and stability of the inhibitor in the fertilizer, which can vary depending on its compatibility with physical mixtures, time, and storage temperature.
The use of stabilized fertilizers may not be the most suitable technology among the strategies to reduce N-NH 3 losses in a brazilian coffee crop system, as they depend on the occurrence of rain to incorporate the N fertilizer.The highest N losses through volatilization for conventional urea occur in the first 7 days after fertilization (Costa et al., 2003;Afshar et al., 2018), which is a critical period of losses (CPL).NBPT typically delays the CPL to increase the chance of incorporation of the N fertilizer into the soil by the precipitation.However, relying solely on weather conditions represents a high risk of inefficiency in N fertilizer incorporation, as it also depends on both the volume and intensity of rainfall (Dawar et al., 2011).In addition, the architecture of the coffee plants creates a barrier that reduces the rainfall under the canopy and the direct incidence on the N fertilizer applied in that location.
Currently, microbiological attributes have been used and disseminated as a tool to ensure ecosystem sustainability and monitoring the health of soil studies.The attributes most used as "microbiological indicators of soil quality" are microbial biomass carbon, microbial biomass nitrogen, microbial basal respiration, and enzymatic activities (Doran and Parkin 1994;Dick et al., 1996).The enzymes most used for this purpose are related to the C, N, P, and S cycles (Gil-Sotres et al., 2005;Moreira and Siqueira, 2006;Bastida et al., 2008).
Our hypothesis was that the technologies for nitrogen fertilizers can help mitigate N losses and improve agronomic efficiency and microbiological attributes in coffee fields.Therefore, the aims of this study were: i) quantify N-NH 3 losses from two conventional sources of N and from urea stabilized with NBPT at two doses of N (150 and 400 kg ha -1 ); ii) determine the effect of nitrogen fertilizer technologies and doses of N (0, 150, 275, 400, and 525 kg ha -1 crop season -1 , divided into three applications) for improvement in the nutrition and yield of coffee plants; iii) update the numbers of N exportation by the beans of recent coffee cultivars and understand the real N demand to avoid nitrate losses by leaching; iv) evaluate the effect of nitrogen fertilization on soil microbiological attributes (microbial biomass carbon, microbial biomass nitrogen, microbial basal respiration, metabolic quotient, urease, b-glucosidase, acid phosphatase and arylsulfatase enzyme activities; and, finally, v) support researchers, commercial farmers and fertilizer industries with a new insight and inedited results of N fertilizer technologies as interesting options for reducing N-NH 3 losses by volatilization and increasing the sustainability of coffee production systems.

Characterization of the experimental area
The experiment was conducted during the 2019/2020 and 2020/ 2021 crop seasons with coffee plants under field conditions.The commercial plantation belongs to the Lagoa Coffee Plantation, owned by the NKG/Fazendas Brasileiras group, located in the municipality of Santo Antonio do Amparo, in the state of Minas Gerais, Brazil (Figure 1).This municipality is located in a traditional coffee-producing region in Brazil within the Campos das Vertentes geographical indication, at 1,100 meters of altitude, latitude 20°53'26.04"S, and longitude 44°52'04.14"W. According to the Köppen classification, the climate is Cwa, with a humid tropical climate with a dry winter and a temperate summer.The mean annual precipitation is approximately 1,493 mm and the mean annual temperature is 19.6°C.
The soil was classified as an Oxisol according to United States Department of Agriculture (USDA) Soil Taxonomy (Soil Survey Staff, 2018).Before the experiment, soil samples were collected at 0 to 20 cm depth to identify the texture as 32%, 8%, and 60% of sand, silt, and clay, respectively, using the Bouyoucos method (Bouyoucos, 1951).
The five years old plantation in the production phase was planted with Coffea arabica L., Catuaı́Vermelho -IAC 99 cultivar.Plants were spaced by 3.40 m between rows and 0.6 m between plants, resulting in 4,902 plants ha -1 .

Fertilizer characterization
The fertilizers used in this study were chosen from the main technologies currently used in Brazilian coffee production.These fertilizers belong to two major technology groups: conventional and stabilized.All fertilizers were photographed using an Olympus microscope, SZ60 Japan model (S.1).For the conventional group, U (46% N) (S.1a) and AN (33% N) (S.1b) were used.For the stabilized group, U NBPT (46% N) (S.1c) was used.
The concentration of NBPT in the fertilizers was analyzed by liquid chromatography (HP1100 Agilent model) with a diode array detector, according to the method described by the European Committee for Standardization (2008).This procedure was done to verify if the concentration of the additive in U NBPT was the same as the manufacturer's specification (530 mg kg -1 of NBPT).However, at the time of application, the NBPT concentration was  Location of the experimental areas in the municipality of Santo Antônio do Amparo, Minas Gerais state, Brazil.
110 mg kg -1 , indicating that there was degradation of the molecule between the manufacturing and commercialization.Some conditions such as storage time and temperature (Cantarella et al., 2018) and contact with phosphate fertilizers and ammonium sulfate that contain free acidity (Sha et al., 2020) accelerate the degradation of this additive and, consequently, reduce the potential for enzyme inhibition.
The treatment of urea with this additive consists of spraying 2.5 to 3 L of a 20% NBPT solution per t -1 of fertilizer.This solution is composed of the NBPT molecule, organic solvent, and dye used to distinguish the stabilized from the conventional fertilizers.

Experimental design
A randomized block design was used in a 3 × 5 factorial scheme, composed of three fertilizers and five doses with four replications, totaling 60 experimental units.Each experimental unit was composed of 16 coffee plants.The area comprising the central 10 plants was used as the useful area for data collection.An entire row on each site was left as border between blocks (Figure 2).
The treatments were composed of U, AN, and U NBPT , applied at doses of 0, 150, 275, 400, and 525 kg ha -1 year -1 of N. The fertilization was split into three equal parts applied 30 to 40 days after the previous fertilization, starting in November, at the beginning of the rainy season.All fertilizers were applied as topdressing and under the projection of the canopy of the coffee plants.

Complementary fertilization management
Liming was done 60 days before the application of the treatments to increase soil base saturation to 60%.Each treatment plot received 1.1 t ha -1 of limestone (90% ENP 16% MgO), 1.5 t ha -1 of gypsum, and 400 kg ha -1 of magnesium oxide (70% MgO).Additional fertilization was done with KCl (60% K 2 O) at 260 kg ha -1 year -1 of K 2 O divided into three applications and triple superphosphate (TSP -46% P 2 O 5 ) at 100 kg P 2 O 5 ha -1 year -1 applied in a single application.Both KCl and TSP were applied in the projection of the coffee tree canopy.
Micronutrients were applied via foliar fertilization along with the plant disease control treatments during the formation and management of the plantation.A commercial product with the following composition was applied at 5 kg ha -1 in 300 L ha -1 of solution: 6.0% of zinc (ZnSO 4 ), 3.0% of boron (H 3 BO 3 ), 2.0% of manganese (MnSO 4 ), 10% of copper (Cu(OH) 2 ), 10% of sulfur (SO 4 ), 1.0% of magnesium (MgSO 4 ) and 10% of potassium (KCl).This fertilization was divided into three applications per year at intervals of 45 days between November and February.The soil was also fertilized with 3 kg ha -1 of B as H 3 BO 4 and 5 kg ha -1 of Zn as ZnSO 4 along with the first application of the maintenance fertilization.

Monitoring of climatic conditions
A climatological station set up near the experimental area was used for daily monitoring of rainfall, temperature (maximum and minimum), and relative humidity throughout the period of conducting the trials in the four agricultural crop seasons.Rain collectors were also installed beneath the coffee tree canopy to assess how much these plants reduce the direct incidence of rain where limestone and fertilizer are applied.

Quantification of the N-NH 3 losses
A 3 × 2 factorial arrangement was adopted involving the same plots with the N sources of experimental design (described above) and two of those doses of N (150 and 400 kg N ha -1 crop season -1 ) were selected.The losses of N through NH 3 volatilization resulting from soil N fertilization were quantified using the semi-open collector, adapted from Lara-Cabezas et al. (1999) method for capturing N-NH 3 .In each experimental plot that received 150 and 400 kg ha -1 of N, three 200 mm diameter bases made of PVC, with 20 cm height, were installed at 5 cm depth in the soil under the projection of the canopies (Figure 3).
PVC chambers similar to the bases but measuring 50 cm height were used to collect the emissions.Specifications are presented in Figure 3. Within each base (0.2 linear meters), the amount of fertilizer corresponding to the dose applied per hectare was added.
Immediately after fertilization, collection chambers were added over the bases in all plots.Within each chamber, two laminated foam disks with 0.02 g cm -3 density and 2 cm thickness, cut to the same diameter as the tube, were inserted.The disks were positioned 30 and 40 cm above the soil.The lower sponge was previously soaked with 80 ml of a solution composed of H 3 PO 4 (60 ml L -1 ) and glycerin (50 ml L -1 ) to retain the volatilized ammonia.The upper sponge served to protect against possible environmental contamination.

FIGURE 3
Illustration of the collectors used for the quantification of ammonia emissions.Sarkis et al. 10.3389/fpls.2023.1291662Frontiers in Plant Science frontiersin.org The lower sponges were collected on the 1 st , 2 nd , 3 rd , 4 th , 5 th , 6 th , 8 th , 10 th , 13 th , 16 th , 20 th , 24 th , 28 th , and 32 nd day after each fertilizer application.After replacing the sponges, the chamber was switched to the other base to reduce the spatial variability of the ammonia emissions.This rotation allowed a greater influence of climatic variations, such as temperature and precipitation.
The solution contained in the sponge was then extracted by filtration through a Buchner funnel and a vacuum pump, along with five consecutive washes with 80 mL of distilled water.From the extract, 20 mL were taken to determine the N content by distillation using the Kjeldahl (1883) method.The N content in the sample was calculated according to Equation 1: in which TN: Total nitrogen concentration in the sample (%); Va: Volume of hydrochloric acid solution used for titration of the sample (mL); Vb: Volume of hydrochloric acid solution used for titration of the blank (mL); F: Correction factor for 0.01 M hydrochloric acid; P1: mass of the sample (g).
The N concentration in the samples collected in the area of the chambers was extrapolated to the percentage of N-NH 3 loss per hectare.To calculate the accumulated losses during the evaluation days, the losses from the 1 st and 2 nd day were added, then the sum of these was added to the 3rd day, and so on until the last day of sampling.

Soil sampling
To quantification of urease activity, soil samples were collected at 2 cm depth from the line of the N fertilizer application on the same dates that the sponges used to capture N-NH 3 were collected.A split-plot design with repeated measures over time was adopted with four replications.The soil, with the same field moisture, was sieved through a 2 mm mesh to remove plant debris, and then stored at 4°C until analysis.
To evaluate the effects of the N sources on microbial C and N, basal respiration and metabolic quotient, b-glucosidase, acid phosphatase, and arylsulfatase, a split-plot design with repeated measures over time was also adopted with four replications.The design consisted of two sampling times in plots with no N (dose 0) and in plots with 400 kg ha -1 of N of AN, U, and U NBPT sources to assess the impact of this dose, often used by farmers, on microbiological attributes in productive and high-investment coffee plantation.Soil samples from the 0-0.10 m layer were collected in the 2021/2022 crop season, both in the projection of the coffee tree canopy and in the fertilization line of each plot.The samples were collected two times: immediately before the start of the N fertilization and after 2/3 (267 kg N ha -1 ) of the N application to evaluate the influence of treatments that have been applied for consecutive years.The samples (2000 g each) were sieved through a 2 mm mesh and stored at 4°C for the following analyses.Other fertilization and management practices were kept constants for all plots, following the same detailed procedure as previously described.

Microbial C and N
Microbial biomass carbon (MBC) and microbial biomass nitrogen (MBN) in the soil were determined in triplicate by the irradiation/extraction method (Islam and Weil, 1998) using 20 g of soil and 50 mL of potassium sulfate (0.5 mol L -1 ).The determination of MBC was carried out through oxidation with potassium dichromate (0.066 mol L -1 ) and subsequent addition of H 2 SO 4 and H 3 PO 4 .The resulting mixture was heated on a hot plate until ebullition.After cooling, the solution was titrated with ammonium ferrous sulfate (0.033 mol L -1 ) and 1% diphenylamine as an indicator (Vance et al., 1987).MBN was determined by the Kjeldahl method.

Basal respiration and metabolic quotient
The soil basal respiration (SBR) was quantified by incubating 20 g of soil in a hermetically closed container with NaOH solution (0.5 mol L -1 ) for three days at 28°C.The reaction was stopped by adding BaCl 2 .2H 2 O (0.5 mol L -1 ) and the solution was titrated with HCl (0.5 mol L -1 ) (Anderson and Domsch, 1993).The metabolic quotient (qCO 2 ) was calculated by the division SBR/MBC (Anderson and Domsch, 1993).

b-glucosidase, acid phosphatase, and arylsulfatase
Samples of 1 g of soil, containing 1 mL of substrates (rnitrophenyl-b-D-glucoside for b-glucosidase, r-nitrophenylphosphate for acid phosphatase and r-nitrophenyl-sulfate for arylsulfatase) were incubated for one hour at 37°C in the presence of toluene and standard solution at specific pH for each activity (pH 6.0 for b-glucosidase, 6.5 for acid phosphatase and 5.8 for arylsulfatase).After the incubation period, the reaction was stopped with the addition of CaCl 2 (0.5 mol L −1 ) and NaOH (0.5 mol L −1 ).The supernatant was filtered and readings were taken in a spectrophotometer at 410 nm by the colorimetric determination of the r-nitrophenol using 4-nitrophenyl b-Dglucopyranoside, potassium 4-nitrophenyl sulfate, and disodium 4-nitrophenyl orthophosphate as substrates, respectively (Dick et al., 1996;Aragão et al., 2020).The values of each activity were expressed in μg PNP g -1 soil h -1 .

Urease activity
The quantification of the urease activity was based on the determination of the N-NH 4 + released after soil incubation with urea solution (Tabatabai and Bremner, 1972).For each treatment, a 5 g soil sample was mixed with 0.2 mL of toluene, 9.0 mL of buffer solution (pH 9.0), and 1.0 mL of a 0.2 mol L -1 urea solution, and then incubated for 2 hours at 25 °C in a temperature and humidity controlled chamber (model TE-371 TECNAL).After incubation, 35 mL of a 100 ppm solution composed of 2 mol L -1 KCl and 2 mol L -1 Ag 2 SO 4 were added to stop the reaction.After shaking, the solution was left to rest for 5 minutes at room temperature, and the volume was then completed to 50 mL with KCl and Ag 2 SO 4 , of which 20 mL were distilled with 0.2 g of magnesium oxide.The distillate was collected in flasks containing boric acid and the indicators methyl red and bromocresol green, and then titrated with a standard 0.005 mol L -1 H 2 SO 4 solution.The control followed the same procedures described above, but the urea was added after the addition of the KCl and Ag 2 SO 4 solution.
2.8 N content in the beans, N content in the leaves, and N exportation by the beans The N content in the beans and leaves were determined by the Kjeldhal method (Kjeldahl, 1883) after sulfuric digestion (Bremner, 1996).Samples of coffee beans were also collected from nearby commercial farms of the same age as the experimental area.These samples belonged to 8 different cultivars, commonly planted in the region.The beans were dried and processed to quantify the bean/ husk ratio and the respective N content, according to the international NIST standard.This procedure allowed the comparison between the proportion of beans/husks and the exported amount of N of the main coffee cultivars, which is a piece of information currently unavailable in the literature.

Yield
The yield of the coffee plants was also evaluated in both crop seasons by manually harvesting the fruits when approximately 70% of them were mature (cherry stage).The harvested volume was measured, and a 5 L sample was sun-dried until reach 11.5% humidity.After drying, the sample was processed, weighed, and subsequently converted into 60 kg bags of processed coffee per hectare.

Calculations and statistical analyses
Data were subjected to a non-linear regression analysis using a logistic model (Equation 2) to evaluate the variable ammonia loss through volatilization.This model is already well-used for estimating plant growth, but it has been recently used to estimate the accumulated N-NH 3 loss (Soares et al., 2012;Silva et al., 2017;Minato et al., 2020).
Where Yi is the i-th observation of the cumulative loss of N-NH 3 in %, being i = 1, 2… n; daai is the i-th day after fertilization; a is the asymptotic value that can be interpreted as the maximum cumulative amount of N-NH 3 loss; b is the abscissa of the inflection point and indicates the day when the maximum loss due to volatilization occurs; k is the value that indicates the precocity index and the higher its value, the lower time needed to reach the maximum loss by volatilization (a); Ei is the random error associated with the i-th observation, which is assumed to be independent and identically distributed according to a normal distribution with mean zero and constant variance, E ~N(0, I s 2 ).
To estimate the maximum daily loss (day of the maximum loss of N-NH 3 ), that is, to determine the inflection point of the curve, the following equation was used: Where k is a relative index used to obtain the maximum daily loss (MDL), and a would be the asymptotic value that can be interpreted as the maximum cumulative loss of N-NH 3 .
After checking for normality (Shapiro-Wilk's test) and homogeneity of variances (Bartlett's test), an analysis of variance was applied to test the influence of the N sources applied on the variables of N-NH 3 loss by volatilization, urease activity, N concentration in the leaves, and yield.The significance of the differences was evaluated at P ≤ 0.05, and after validation of the statistical model, means were grouped by the Scott-Knott's algorithm using R software 3.3.1 (R Development Core Team, 2018).

Weather conditions
The accumulated precipitation on the 1 st , 2 nd , and 3 rd applications of the N fertilizer in the 2019/2020 crop season were 289, 457, and 336 mm, respectively, totalling 1082 mm.In the first seven days after each fertilization, the precipitation corresponded to 0 mm (0%), 63 mm (14%), and 159 mm (47%).
In the 2020/2021 crop season, the accumulated precipitation on the 1 st , 2 nd , and 3 rd fertilization were 291, 280, and 336 mm, respectively, totalling 907 mm.In the first seven days after each fertilization, the precipitation corresponded to 14 mm (15%), 28 mm (10%), and 0.2 mm (0.05%).The precipitation influenced the processes related to N-NH 3 losses and we observed that lower rainfall in the first seven days after fertilization resulted in lower N-NH 3 losses.
The area under the coffee canopy, where liming and fertilizers are applied, received 30 to 40% less precipitation compared to the inter-row area, depending on the rainfall volume and intensity, plant architecture, leaf area, and development stage.
The relative air humidity was above 70% for most of the period following the N fertilization in both crop seasons.The average temperature during the same period was 22 °C, with a minimum of 18.3 °C and a maximum of 36.6 °C.

Daily and cumulative losses of N-NH 3
Daily and cumulative losses of N-NH 3 from the three fertilizer applications in each crop year were influenced (P ≤ 0.05) by the N fertilizer technologies.The volatilization of N-NH 3 was similar for the doses of 150 and 400 kg ha -1 of N. The application of U resulted in losses of 22.0% and 22.8% for these respective doses, considering the average of the two evaluated crop seasons.This means that 66 and 182 kg ha -1 of N were lost, respectively, at the end of six fertilizations with conventional urea (Table 1).U NBPT reduced N-NH 3 losses compared to U. Average losses of 16.9% and 19.8% were observed for doses of 150 and 400 kg ha -1 of N supplied by this stabilized source with a urease inhibitor additive.This reduction in volatilization prevented the loss of 15.9 and 24.3 kg ha -1 of N, respectively, compared to U (Table 2).As for AN, N-NH 3 losses did not exceed 1% of the applied dose, regardless of the weather conditions at the time of the fertilization.This reduction in volatilization was 90.9 and 95.8% lower, respectively, preventing losses of 60 and 174.8 kg ha -1 of N compared to conventional urea (Table 1).
In the 2019/2020 crop season, the maximum daily N-NH 3 volatilization for U occurred 1.4 (7.36 kg ha -1 of N) and 1.88 (5.11 kg ha -1 of N) days after the application of 150 and 400 kg ha -1 of N. For AN, the maximum daily losses of 0.05 and 0.09 kg ha -1 of N occurred 7.31 and 5.76 days after the fertilization, respectively.For urea stabilized with NBPT, the maximum daily losses of 4.33 and 2.38 kg ha -1 of N occurred at 1.9 and 2.04 days after fertilization, respectively (Table 3).
In the 2020/2021 crop season, the maximum daily loss for U occurred 1.77 and 3.11 days after the N fertilization, with 9.66 and 4.07 kg ha -1 of N, respectively.For AN, the maximum daily loss occurred 9.49 and 12.90 days after fertilization, with similar volatilization to the first season (less than 0.5 kg ha -1 of N).U NBPT presented maximum losses at 3.48 and 2.82 days, with the volatilization of 4.55 and 6.37 kg ha -1 of N, respectively (Table 2).

Urease activity
Urease activity was higher on days of maximum NH 3 volatilization.In general, the following sequence for the treatments was observed: U > U NBPT > AN, despite some variations that are common for these factors (Figure 4).U NBPT reduced the urease activity, proving that the additive associated with the fertilizer prolonged the urea hydrolysis due to the enzymatic inhibition that resulted in a decrease in N-NH 3 losses presented in item 3.2.

Microbiological attributes
There was an effect of the N sources (NS), soil sampling time (ST), and their interaction (NS × ST) on the MBN and arylsulfatase activity.The N sources also influenced the MBC and the metabolic quotient (qCO 2 ).
In general, the soil showed high biological quality.The average content of N in the soil microbial biomass was higher after the N fertilization (average MBN = 644 mg N g -1 of dry soil), and the U NBPT stood out (1681 mg N g -1 of dry soil) (Figure 5A).The MBN of this treatment before the fertilization (45 mg N g -1 of dry soil) was lower than the control without N fertilization (216 mg N g -1 of dry soil) and in the soil fertilized with conventional urea (354 mg N g -1 of dry soil).Soil microbial C content varied significantly with the N sources only before the fertilization (48 to 430 mg C g -1 of dry soil), where it decreased in the following order: U NBPT = control > ammonium nitrate > conventional urea (Figure 5B).The sampling time did not interfere with the MBC in the control and urea + NBPT treatments.
The average SBR was 12 mg C-CO 2 kg -1 dry soil h -1 .The results of basal respiration of the N fertilizers remained equivalent to the control treatment in both sampling times (Figure 5C).Conversely, the soil fertilized with conventional urea had the highest qCO 2 before the fertilization (Figure 5D).This result is directly related to the low microbial carbon identified in the initial sampling (48 mg C g -1 dry soil).The high qCO 2 with conventional urea (321 mg h -1 C-CO 2 mg -1 C-mic) represents about seven times the average of the other treatments for the same sampling time.However, this effect was not observed in the qCO 2 after fertilization, where the result with conventional urea was statistically similar to the other treatments (Figure 5D), thus indicating the variability of these parameters.
The activities of b-glucosidase and acid phosphatase were not influenced by the addition of the N fertilizers in either of the soil sampling.High levels of p-nitrophenol (PNP), which reflect the enzymatic activity, were observed in both sampling times.For bglucosidase (Figure 6A), PNP activity ranged from 3780 to 4970 μg PNP g -1 of dry soil before the fertilization and from 2260 to 4550 μg PNP g -1 of dry soil after the fertilization.For acid phosphatase (Figure 6B), PNP activity ranged from 660 to 1538 μg PNP g -1 of dry soil before fertilization and from 620 to 1405 μg PNP g -1 of dry soil after fertilization.
The arylsulfatase activity before fertilization was higher in the soil with NBPT, while the other N sources did not differ significantly in the first sampling time (Figure 6C).After fertilization, the control results were superior, especially compared to the conventional urea, which represented only 18% of the total activity of the treatment without N fertilization.There was a significant decrease in arylsulfatase activity between the samplings, with reductions of 45%, 69%, and 77% in the levels of p-nitrophenol (PNP) immediately after fertilization with ammonium nitrate, urea + NBPT, and conventional urea, respectively.

N content in the beans and leaves
The total N content in the beans of the 2019/2020 and 2020/ 2021 crop seasons were not influenced (P ≤ 0.05) by the isolated effect of the N doses and sources, and there was no significance (P ≤ 0.05) for the interaction between these factors (Table 4).The average N content was 2.6% in the beans and 1.6% in the husk, Soil urease activity after N fertilizer application with split fertilization (0 and 400 kg ha -1 of N) in coffee plantation.
while the average bean/husks ratio is 1:1, despite the cultivar or productive load.
The average N content in the leaves (30 mg kg -1 ) of the samples collected during the green cherry stage was not influenced (P ≤ 0.05) by the isolated effect of the N doses and sources.There was also no significance (P ≤ 0.05) for the interaction between these factors (Table 5), which is explained by the high levels of N present in the soil.
The extraction and export of N by coffee plants depends on the crop's productivity and vegetation intensity, which vary in accordance with the biennial nature of the crop and annual soil and climate conditions.

Yield
Treatments did not lead to different yields (P ≤ 0.05) in the 2019/2020 crop season.For the 2020/2021 crop season, there was a significant effect (P ≤ 0.05) of the N sources within the dose of 525 kg ha -1 , with AN being higher than U and U NBPT (Table 6).
In both seasons, no relationship was observed between the amount of N lost by volatilization and yield, meaning that yield was not affected even in the treatments that showed higher losses by volatilization (Table 6).

Weather conditions, N-NH 3 volatilization and urease activity
In this study, the climatic conditions, especially precipitation and temperature, directly influenced the losses of N-NH 3 by volatilization.In both coffee seasons, most of the N-NH 3 losses occurred during the first seven days after the split application of the N fertilizers.The precipitation during these first days was essential to incorporate the fertilizer into the soil and reduce the N-NH 3 .
Rainfall, depending on the volume and intensity, enables the incorporation of the N fertilizer into the soil (Dawar et al., 2011).Higher rainfall volume and intensity are usually associated with better fertilizer incorporation and, consequently, lower losses by volatilization.The lower precipitation volume immediately after the application of the fertilizers on the soil hindered the incorporation of the N and favored N-NH 3 losses.Thus, our results indicated that relying solely on environmental conditions for fertilizer incorporation increases the risk of N losses and, consequently, lead to the low efficiency of the fertilization.
Moreover, the architecture of the coffee plant obstructs the direct incidence of rainfall that would incorporate the N fertilizer applied under the canopy.Additionally, the presence of plant residues on the soil also acts as a barrier to fertilizer incorporation and creates a favorable environment for volatilization due to the high concentration of urease.The low concentration of NBPT demonstrated by chromatography analysis, as mentioned, may often not be sufficient to inhibit urease activity, resulting in daily N-NH 3 losses similar to those observed for conventional urea.
Ammonium nitrate is a widely used N fertilizer in Brazilian coffee cultivation.In this study, the use of this fertilizer reduced over 90% of the N-NH 3 losses in both crop seasons.The insignificant loss of N-NH 3 for this N source is related to its acid to neutral reaction in the soil, especially at pH< 7. Another positive aspect is that these fertilizers do not depend on weather conditions at the time of application (Freitas et al., 2022;de Souza et al., 2023).Therefore, ammonium nitrate can be a strategic option for N fertilization in coffee cultivation systems.Conversely, amidic N is reported as the most susceptible to N-NH 3 losses and represents an average of 30% of the N applied in coffee cultivation areas in Brazil (Chagas et al., 2016;Dominghetti et al., 2016;Silva et al., 2017).This means that for every three fertilizer applications, approximately one is lost due to volatilization.For an average annual fertilization of 845 kg of urea per hectare (380 kg of N ha -1 ), 254 kg of urea per hectare (114 kg of N ha -1 ) are lost.This currently corresponds to a financial loss of US$ 228 per hectare (FL= average annual loss of N-NH 3 x average value per kg N-fertilizer), considering US$897 as the recent average value per ton of fertilizer.For ammonium nitrate, these losses rarely exceed 1%, an amount that is probably emitted by the soil, which justifies the greater retention and storage of N in soil fertilized in the correct dose, time, and location.

Microbial and enzymatic activity
The soil microbiological attributes indicated high biological quality for this adequately fertilized scenario of productive coffee farming.The results demonstrate that soil sampling over time (before or after fertilization) can identify different levels of microbial activity, highlighting spatial and temporal variability for soil quality indicators.These indicators are not always related to the fertilization carried out in that crop season, but to the fertility management adopted over several years that favor microbial resilience and sustainability of the productive system.
The abundance of this enzyme in the soil is an indicator of a high rate of organic sulfur mineralization and high fungal biomass, which is the main producer of soil sulfate esters (Dick et al., 1996).Some factors such as soil texture, organic matter, moisture, and oxygen affect the survival of microorganisms (Zantua and Bremner, 1977;Burns, 1982;Flores-Renterıá et al., 2020) and, consequently, the level of microbial and enzymatic activity.Thus, soil biological quality indicators have little contribution to decision-making regarding short-term agronomic management of coffee crops, and the results reflect the management adopted for several years, necessary for significant changes to occur in C stocks (Sarkis et al., 2021), soil profile building for adequate plant development, organic residue input, and physical, chemical, and biological balance of the soil.

N content in the beans, N content in the leaves, and N exportation by the beans
The N content in coffee beans is a physiological characteristic of the plants themselves, independent of the cultivar (Table 5).The green cherry stage is the time of highest plant demand (Prieto Martinez et al., and this plants redistributes the stored N to the fruits, maintaining these levels almost constant, regardless of the productive load.In addition to the N exported by the processed coffee and the absorption necessary for plant growth (Table 4), a large amount of N is removed from the system by the husk of the harvested fruits, which represents approximately 50% of the harvested fruit weight (Tables 4, 5).Is emphasized that the exceeding doses of N applied in coffee plantations results in a N mineral movement from the surface to the subsurface (N loss by leaching) despite the source of the fertilizer (Sarkis et al., 2023).Cannavo et al. (2013) also observed that the addition of large amounts of N fertilizer every year in this system led to a soil N saturation and a high potential of N losses through NO 3 -.

Yield
In the experiment, the reduction in ammonia volatilization did not influence crop yield.Other authors reported that N volatilization, in the form of NH 3 , does not limit production (Chagas et al., 2016).This may be explained by the amount of N applied as fertilizer being more than sufficient for plant growth.The biennial cycle of coffee, with one year being more productive than the following, can make difficult the observation of the effects of treatments on yield.Furthermore, the N stock in the soil organic matter should be taken into account, as it can compensate for the N lost through volatilization in conventional fertilizers, providing a good short-term productivity for the crop.However, over time, there may be a reduction in the N stocks due to lower N retention in the soil.The lower N-NH 3 losses from the use of AN can increase the average mineral N content by 50% when compared to urea and U NBPT (Sarkis et al., 2021).

Conclusions
Our study showed that N fertilizer technologies are interesting options for reducing N-NH 3 losses by volatilization, increasing N retention in the soil, and the sustainability of coffee production systems.U NBPT inhibited urease activity in the soil and reduced N-NH 3 losses compared to conventional urea, but it is not the most strategic option in the context of coffee farming because it depends on the volume and intensity of rain to incorporate the fertilizer after the fertilization.Regardless of the weather conditions, AN lost less than 1% of N by volatilization, being the safest option for N fertilization efficiency in coffee production.This characteristic is especially advantageous considering the specific microclimate in the projection of plant canopies that hinders the incorporation of the fertilizers by the direct incidence of rainfall.Thus, AN is the source that contributes the most to increasing soil N stocks and supplying the demands of the vegetative and reproductive stages of the coffee plants.Soil microbial and enzymatic activity indicated high biological quality for this appropriately fertilized productive coffee farming.NBPT, N-(n butyl) thiophosphoric triamide; Ept: N export per bag of processed coffee.In each crop season, 150 or 400 kg N ha −1 per year were split into three equal applications for conventional and stabilized N fertilizers, totaling 300 or 800 kg N ha -1 for both crop seasons, respectively.*Means did not differ by the Scott-Knott's test (P ≤ 0.05).Sarkis et al. 10.3389/fpls.2023.1291662Frontiers in Plant Science frontiersin.org Although indicators of soil biological quality have been widely discussed in recent scientific publications, our indicated that these guidelines have done little to aid decision-making regarding agronomic management of coffee crop, reflecting the management practices adopted for several years that have resulted in the accumulation of organic carbon, development of soil profile for adequate plant growth, organic residue input and maintenance of physical, chemical, and biological balance of the soil.

FIGURE 2
FIGURE 2Illustration of the experimental design.

a:
Asymptotic value (percentage of the estimated maximum accumulated loss); b: Day when the maximum ammonia loss occurs; k: relative index; MDL (maximum daily loss of ammonia); and NBPT, N-(n butyl) thiophosphoric triamide.

FIGURE 4
FIGURE 4 FIGURE 5Microbial biomass nitrogen (A) and carbon (B); soil basal respiration (C); and metabolic quotient (D) as a function of the N sources (ammonium nitrate, conventional urea, urea + N-(n-butyl) thiophosphoric triamide (NBPT), and control without mineral N fertilization for six years), for the doses 0 and 400 kg ha -1 of N, and soil sampling time (samples 1 and 2).Means followed by the same uppercase letter between N sources and lowercase letter between soil sampling time do not differ statistically by Scott-Knott's and F tests (P ≤ 0.05), respectively.Standard error bar; mean of four replicates.
FIGURE 6Enzymatic activity of b-glucosidase (A), acid phosphatase (B), and arylsulfatase (C) as a function of N sources (ammonium nitrate, conventional urea, urea + N-(n-butyl) thiophosphoric triamide (NBPT), and control with no mineral N-fertilization for six years), for the doses 0 and 400 kg ha -1 of N, and soil sampling time (sampling 1 and 2).Means followed by the same uppercase letter among nitrogen sources and lowercase letter among soil sampling time do not differ statistically by Scott-Knott's and F tests (P ≤ 0.05), respectively.Standard error bars; mean of four replicates.PNP = r-nitrophenol.

TABLE 1 Mean
Cumulative losses of ammonia (% of the applied N), for conventional and stabilized N fertilizers, in three fertilizations in the coffee plantation, during the 2019/2020 and 2020/2021 crop seasons.In each crop season, 150 or 400 kg N ha −1 per year were split into three equal applications for conventional and stabilized N fertilizers, totaling 300 or 800 kg N ha −1 for both crop seasons, respectively.*Means followed by the same lowercase letter in the column do not differ by the Scott-Knott's test (P ≤ 0.05).Mean of the six fertilization sources performed between November and February during both seasons.

TABLE 2
Adjusted regression parameters for the accumulated and maximum daily losses of N-NH 3 from conventional and stabilized N fertilizers in the 2020/2021 crop season.
a: Asymptotic value (percentage of the estimated maximum accumulated loss); b: Day when the maximum ammonia loss occurs; k: relative index; MDL (maximum daily loss of ammonia); and NBPT, N-(n butyl) thiophosphoric triamide.

TABLE 3
Adjusted regression parameters for the cumulative and maximum daily losses of N-NH 3 from conventional and stabilized N fertilizers in the 2019/2020 crop season.

TABLE 4 N
content in the beans (%) and N export per bag of processed coffee (kg N bag -1 ) for conventional and stabilized N fertilizers in a coffee plantation, during the 2019/2020 and 2020/2021 crop seasons.

TABLE 6
Coffee bean yield (bags ha -1 ) according to the interaction between the N sources and N doses during the 2019/2020 and 2020/2021 crop seasons., Ammonium nitrate; U NBPT , urea + N-(n-butyl) thiophosphoric triamide (NBPT).Uppercase letters cluster sources (with four replications each of each one of the doses) according to the Scott-Knott's test. AN

TABLE 5
Percentage of beans and husks, and N content in the beans of different coffee cultivars at full productive potential, harvested during the 2020/2021 crop seasons.
NIST, International standard for validation of foliar analysis.